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Abstract. In this work we study the high pressure behaviour of liquid and glassy 
Ge02 by means of molecular dynamics simulations. The interaction potential, 
which includes dipole polarization effects, was parameterized from first-principles 
calculations. Our simulations reproduce the most recent experimental structural 
data very well. The character of the pressure induced structural transition in the 
glassy system has been a matter of controversy. We show that our simulations and 
the experimental data are consistent with a smooth transition from a tetrahedral to 
octahedral network with a significant number of penta-coordinated germanium ions 
appearing over an extended pressure range. Finally, the study of high-pressure, liquid 
germania confirms that this material presents an anomalous behaviour of the diffusivity 
as observed in analogous systems such as silica and water. The importance of penta- 
coordinated germanium ions for such behaviour is stressed. 
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1. Introduction 

Germania (Ge02), along with silica and beryllium fluoride, is a "strong" glass-former [1] 
characterized by a tetrahedral network structure at ambient conditions in the amorphous 
phase p]. Its structure is based on corner-sharing Ge(Oi/2)4 tetrahedra (the 1/2 index 
means that each ion is shared by two Ge^'^), with a Ge-0 average distance of 
1.73 A and a mean inter-tetrahedral Ge-O-Ge angle of 132° [3lll]. There is considerable 
interest in the behaviour of such tetrahedrally coordinated glasses under pressure. In 
the associated crystalline materials high-pressure transitions are observed between four- 
and six-coordinate structures. Pressure-induced structural changes in the amorphous 
phases have been linked to anomalous behaviour in the elastic, viscous and thermal 
properties and to the phenomenon of polyamorphism [5] . In amorphous silica itself, the 
major structural transition occurs at relatively high pressures (~25 GPa) where direct 
structural studies remain difficult [6J. However, because of the larger cation / anion 
radius ratio in germania relative to silica, the transition occurs in a pressure domain 
now accessible to structural studies. Despite this, a clear picture of the nature of the 
pressure-induced structural changes in germania has only recently started to emerge. We 
have therefore performed computer simulations in an attempt to clarify the relationship 
of the information emerging from the different experiments. 

High-pressure structural studies have been made using EXAFS spectroscopy [71 El 
[9l [To] and X-ray and neutron diffraction experiments [U [TTl [12]. All the experimental 
studies agreed on the existence of a structural change in the pressure range 3-15 GPa, 
associated with an increase in the Ge-0 separation which is broadly consistent with an 
increase in coordination number. Vaccari et al. [8], for example, showed, from EXAFS 
studies, that this distance switched from 1.74 A (at GPa) to 1.82 A (at 13 GPa). On 
the other hand, the literature has long been very contradictory about the character of 
the change in coordination number. The first study to address this issue is that of Itie 
et al. [7] who performed x-ray absorption measurements up to 23.2 GPa and reported 
that the Ge coordination changes from fourfold to sixfold at pressure between 7 and 9 
GPa. Vaccari et al., on the other hand, proposed a progressive shift in the coordination 
number, and found no evidence of a fully six-coordinate structure, even at the highest 
pressure of their study (13 ~ GPa). Very recent XAFS and EXAFS studies P [ID] 
extended the studied pressure range to 53 and 44 GPa respectively and it was postulated 
that a complete 6-fold coordination of the Ge ions is only reached at pressures as high as 
25-30 GPa. The ability of EXAFS to provide very accurate first-neighbour distances is 
well established, but so is its limitation for the determination of coordination numbers 
in amorphous materials [131 [lH US]- Only an average coordination number can be 
extracted, with an important error bar which depends crucially on assumptions about 
the shape of the radial distribution functions at larger separations, and to estimate 
proportions of different coordination polyhedra in a material from EXAFS alone is not 
normally considered reliable. 

In principle neutron diffraction, when employing the isotopic substitution method. 
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gives access to all the partial structure factors, and hence to the corresponding partial 
radial distribution functions (RDF) which contain all the structural information. In the 
case of Ge02, this programme has been fulfilled only at ambient pressure by Salmon 
et al. [16]: at elevated pressures only total X-ray and neutron diffraction patterns are 
available. A first attempt to study the high pressure system by these means was reported 
by Guthrie et al [11]: unfortunately, the data were too noisy to extract good RDFs. 
However, a very sharp transition from a tetrahedral to an octahedral structure was 
proposed, with a small range of pressure in which the system is characterized by five- 
fold coordination around Ge^+ ions. Very recently, Drewitt et al. have performed 
new neutron diffraction measurements up to 8.6 GPa and obtained data of very high 
quality [12]. Their results show a gradual change of the intermediate range order with 
increasing density as manifested by an increase in position and reduction in height of 
the first sharp diffraction peak in the total structure factor. From their data there is 
no evidence in support of an abrupt transformation of the network structure over the 
investigated pressure range, in agreement with the most recent EXAFS experiments 

Computer simulations could help in resolving these differences since they can 
provide a detailed picture of the structure of a system. Unfortunately, to date, all 
the molecular dynamics (MD) calculations involving pressurized germania have been 
performed with pair potentials of limited accuracy [T71|T8|[T9|[2nl[2T|[22|[28|[2il[25]. 
For example, for the most commonly used potential, developed by Oeffner and Elliot [26] , 
two different parameter sets were proposed: an original one, which was fitted from an 
ab initio energy surface, and a so-called rescaled one, which was developed from the 
previous one in order to give a better reproduction of the vibrational properties. For 
this reason, some apparent inconsistencies have arisen in the literature because different 
classical potentials were being used [T71 [231 1211 [28]. Finally, since the OE potential is 
not able to reproduce with a single set of parameters all the ambient pressure properties 
of amorphous germania, it is also reasonable to suspect that it might fail when studying 
the pressure behaviour of this material. This might explain the disagreement with the 
experimental trends found in some papers in the literature. For instance, in ref. [25] , 
the Ge-0 bond distance is found to remain constant in the 3-25 GPa pressure interval, 
whilst in ref. [20] the Ge-0 coordination number is found to increase continuously from 
the lowest pressure (IGPa), both trends being in strong contrast with the experimental 
evidence [71[8l[II]. 

Recently, first-principles molecular dynamics (FPMD) studies of glassy Ge02 have 
also been reported [23l |29l |30l [31]. In principle, the amount of empirical information 
needed to set up a first-principles calculation is minimal and it would normally be the 
method of choice to study the physico-chemical properties of condensed phase systems. 
However, it is computationally very expensive compared with classical molecular dynam- 
ics, which is a major drawback when dealing with glassy systems, where the structural 
relaxation times are necessarily long. These studies have therefore focused on the study 
of the structural and vibrational properties of the glass; only one of them included re- 
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suits for high pressure systems [21] • 

We recently introduced a new interaction potential, which includes many-body 
polarization effects, and succeeded in reproducing the structural and vibrational 
properties of glassy germania at ambient pressure as well as the dynamical properties 
of liquid germania in the 3600-5000 K temperature range [28j, with a single set of 
parameters. This potential was parameterized with reference to extensive Density 
Functional Theory (DFT) electronic structure calculations on disordered condensed- 
phase configurations. Also, similarly constructed potentials have been shown to provide 
an excellent, transferable description of a number of complex oxides [321 [33]. The aim 
of this paper is to show how the use of this potential together with a direct comparison 
with the available experimental data, can improve our understanding of the behavior 
of glassy germania at high pressure. In the end, we will emphasize that the various 
experimental studies are in accord, even if the initial interpretations that were made of 
them did not. 

Once the question of structure is resolved, it is of interest to study the influence 
of structural changes on the physical properties of the system. A well-known effect is 
the existence of dynamical anomaly in silica-based systems when the pressure increases 
[5]. Unlike most other systems which, under compression, tend to show a decreasing 
mobility of the species, in liquid silica the diffusion coefficients increase until they reach 
a maximum for a given pressure, and then decrease. This anomaly was linked to the 
existence of pentacoordinated species [M]. The only evidence for such a behaviour 
in liquid germania was provided by Sharma et al. who showed a decrease of the 
viscosity with pressure, though the pressure range of this study was very limited (0- 
1 GPa) [35j. There are also some simulations studies on high-pressure liquid germania 
but, unfortunately, the results are contradictory. Hoang et al. [19j do indeed observe 
an anomalous behaviour of the diffusivity in Ge02. However, Hung and co-workers 
[21] . who also studied this system by means of classical molecular dynamics, found 
no evidence of such a phenomenon, but again both of these studies involved classical 
pairwise additive interaction potentials of limited accuracy. To address this issue we 
have therefore studied the structural and dynamical properties of liquid germania, at a 
temperature of 4000 K and over a wide range of pressures. 

2. Simulation details 

The interaction potential used in this study has already been described elsewhere [28j . 
The ionic species carry their valence charges (Ge^"^ and 0^~), and the polarization 
effects that result from the induction of dipoles on the oxide ions are accounted for. The 
parameters of the interaction potential for germania were obtained from the application 
of a force- and dipole-matching procedure aimed at reproducing a large set of first- 
principles (DFT) reference data [36] on the condensed phase. 

To obtain structural properties in the glassy state, we performed MD simulations 
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in the NPT ensemble using the method introduced by Martyna et al. [3Z]. Details on 
the generation of a compressed glassy state will be reported below. When dynamical 
properties were computed, the system was simulated in the NVT ensemble using the 
Nose-Hoover chain thermostat method [38]. A relaxation time of 10 ps was used for 
the thermostats and the barostats. The simulations on glassy Ge02 (T = 300 K) 
were performed on a simulation cell containing 432 atoms whereas the NVT simulation 
cells on liquid germania (T = 4000 K) contained 600 atoms. We used a time step of 
1 fs to integrate the equations of motion, and minimization of the polarization energy 
was carried out at each time step using a conjugate gradient method. All simulations 
were equilibrated for at least 0.5 ns and then a subsequent run of 0.5 ns was made to 
accumulate enough statistics. 

3. Compressing a glass: some limitations 

Generating a glassy state by means of computer simulations is a challenging operation. 
Indeed the relatively short time-scales available (a few ns in our case) force us to use 
unrealistically fast cooling rates. This implies that sometimes the glass we obtained 
is further away from equilibrium than the experimental one. In fact, even experimen- 
tally, glasses prepared in different ways show different properties, implying that even 
on a time-scale of hours (i.e. 12 orders of magnitude higher than what we can afford 
with MD simulations) there are still some relaxation effects. This problem becomes 
even more significant when one tries to compress a glass. To this end, Scandolo et 
al. described two different compression methodologies for Si02 [39]. The first one is a 
cold compression route consisting of increasing pressure slowly at ambient temperature 
while the second one, which they called a 'quench-from-the-melt' procedure, consists of 
obtaining the compressed glass by a slow cooling of a compressed sample from a tem- 
perature where atomic diffusion is observable on the time scale of the MD simulation. 
By comparing their results with the available experimental data on compressed glass, 
they showed that the samples obtained with the two procedures are representative of 
the experimental in situ compressed (as, for example, in refs. [HI [9]) and densified forms 
of glass [10], respectively. They found that structural differences between annealed and 
cold-compressed forms were most noticeable at distances of 3.5-4 A. 

In the case of Ge02 there has been little effort so far in this direction. Most of the 
previous work on this system used the cold compression route but a direct quantitative 
comparison with experimental data was never attempted. In this work we tried both 
routes used by Scandolo et al. and compared the results with the experimental density 
vs pressure data. This is shown in figure [H It can be appreciated that neither of the 
two routes reproduces quantitatively the experimental data. In both cases the simu- 
lated glass seems to be less responsive to compression, which we attribute to the limited 
simulation times (see below). It seems, however, that the quench- from-a-melt procedure 
gives a closer agreement with the experimental data than the cold-compressed route. 
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Figure 1. The mass density p for Ge02 glass as a function of pressure for cold 
compression and quench-from-a-melt procedures. Experimental data are from the in- 
situ compression studies of Hong et al. [JT] ((red squares)), Smith et al. [H] (empty 
circles) and Tsiok et al. jj^ (blue broken curve). 



once again indicating that a certain degree of diffusion helps the glass relaxation. It has 
to be remembered, however, that Scandolo et al. showed how the quench-from-a-melt 
route is more representative of the densified forms of the glass, so that a comparison 
with the experimental data on in situ compressed glass, such as in refs. [Ill [IB] might 
not work as well. 

It should be evident from the above discussion that we are unable to reproduce 
quantitatively the equation of state of glassy germania. This should not necessarily 
be taken as a limitation of our potential but as a consequence of the relatively short 
time-scales available in computer simulations. For this reason, in the remainder of this 
paper, we will simulate cold-compressed Ge02 and report our data as a function of 
density, instead of pressure. When a comparison with experimental data is required, 
we will convert the experimental pressures into densities by using the data in figure [1] 
from ref. [1^ l¥T] . In the case of liquid germania, since the the system is very diffusive 
and reaches equilibrium in a few picoseonds, we do not anticipate any problem with the 
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equation of state and we will therefore report our data as a function of pressure. 
4. Glassy germania 

In the top panel of figure [2] we report the average Ge-0 distanc^ as a function of 
increasing density. It is compared to the values obtained by Vaccari et al. ^ and Drewitt 
et al. [12] via EXAFS and neutron diffraction measurements. The agreement is quite 
good, considering the scatter of the experimental data and the intrinsic error associated 
with the pressure-to-density conversion. The simulations, like the experiments, show 
first a Ge-0 bond distance which is fairly constant in the 3.66 - 4.5 g/cm^ density 
range (corresponding to an experimental pressure range of 0-5 GPa) and then a gradual 
increase in the Ge-0 bond distance, starting at a density of 4.5 g/cm^. 

We can interpret that data in terms of compression mechanisms. When a 
tetrahedral system like Ge02 is compressed, three types of structural change could occur. 
Firstly, the system could keep exactly the same structure, but adapt by allowing the 
bond lengths to decrease: obviously this is not the case here. Secondly, the tetrahedral 
units could be kept unmodified, in which case the decrease of accessible volume would 
imply a rearrangement of the network structure, progressively minimizing the size of 
the voids [HI HO] and resulting in changes in the first-sharp diffraction peak. Finally 
the tetrahedral structure could be lost, and higher-coordinated structural units could be 
formed. Such an evolution may induce an increase of the Ge-0 bond length in order to 
allow more than 4 oxide ions into the germanium first coordination shell. The analysis of 
bond length data gives evidence of a succession of the two latter mechanisms, at low and 
high pressures. This was also confirmed by an analysis of the bond angle distributions 
(not shown). 

To verify this, we have determined the germanium coordination number. The 
cut-off used for the coordination number calculation was 2.38 A which corresponds 
to the first minimum of the Ge-0 partial RDF. Again we compare these values with 
the ones obtained by Vaccari et al. [8] and Drewitt et al. [12]. Good agreement is 
found, especially with the neutron data. The simulated coordination number remains 
constant, at a value of 4, until a density of 4.5 g/cm^. In this first regime, germania 
is therefore keeping its tetrahedral structure under compression, but the tetrahedra 
reorient themselves. Then, in a second regime, more highly coordinated germanium 
ions begin to be observed. In the inset of the figure, the percentage of Ge04, GeOs and 
GeOe units are given as a function of density. Once again, at a density of 4.5 g/cm^, 
some Ge^"*" ions start to accept 5 0^~ ions in their vicinity. Then, at p ~ 5.0 g/cm^ 
(corresponding to an experimental pressure of about 7-8 GPa), GeOe units begin to be 
formed. At the highest density studied experimentally, p ~ 5.85 g/civ? (corresponding 
to approximately 15 GPa), a non-negligible proportion of Ge04 units is still observed. 

t This distance was defined as doeO = '', where efr) is the Ge-0 radial distribution function 

Jo ''' S^Vjdr °^ ' 

and icut is the position of the first minimum in the g(r) 
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Figure 2. Top panel: comparison between the simulated (solid line) and experimental 
Ge-0 bond distance from Vaccari [8] (black empty dots) and Drewitt [12] (black filled 
squares) as a function of increasing density. Bottom panel: comparison between the 
simulated (solid line) and experimental, (black empty dots) Vaccari [8] and (black filled 
squares) Drewitt [12], first shell Ge-0 coordination number. The inset in the bottom 
panel shows the density dependence of the percentage concentration of Ge04, GeOs, 
GeOe units. 
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We find no evidence of a state in wliicli tliere are only 5-fold coordinated Ge ions in 
the studied pressure range. These fivefold-coordinated units play an important role in 
liquid germania as we shall show in the next section. Our results can also be compared 
to the recent FPMD simulations on the same system [31]. In the latter, a much sharper 
change in the first-neighbour distances is observed upon compression: At a density of 
5.4 g.cm~^, the Ge-0 distance already reaches a value of 1.89 A, showing important 
differences with experimental data. This is very likely due to the small size of the 
samples (108 atoms) and to the higher quenching rates that have to be involved in 
FPMD simulations. 

At first sight, our results differ significantly from the conclusions of Guthrie et 
al. from their combined X-ray and neutron diffraction study of glassy germania [11] . 
Qualitatively, they observed an analogous mechanism, but with a much sharper increase 
in the coordination number. They proposed a complete transition to a GeOs units based 
structure for pressures ranging from 6 to 10 GPa, and a final completely octahedral 
structure at a pressure of 15 GPa. These values were extracted from a Fourier transform 
of the X-ray diffraction total structure factor (their neutron data being too noisy to allow 
an attempt to extract partials). To show that the discrepancy observed is due to the 
analysis of the X-ray data, rather than a deficiency of the representation of the structure 
in the simulations, we have computed the X-ray structure factor. This was obtained 
using the following relationship 

= CGefGe{qf[SGe-Ge{q) " 1] + '^CQeCo fceiq) fo{q)[SGe-o{q)] + 

clfo{q?[So-o{q) - 1] (1) 

where and fa{q) represent respectively the atomic fraction and X-ray form factor of 
element a. Sa/siq) is the partial structure factor which can be obtained via a Fourier 
transform of the partial radial distribution functions: 

S{q)ap = 6^p + p r 47rr2^^^[(?,^(r) - l]dr. (2) 

Jo q^ 

The g-dependent X-ray form factors are calculated from 

4 

fM) = '^aa,iexp [-ba^i{q/AnY] + d^. (3) 

i=l 

The parameters we used for the calculation of these form factors factors are reported 
in table [H note they are consistent with an ionic representation of the distribution of 
electrons {i.e. Ge^"*" and O^^). Figure |3]has been organised in a similar way to Fig. 1 
in ref. [TT] in order to facilitate a comparison; in particular, S{q) has been normalized 
by: 

2 

{Y.^.Uq)f- (4) 

a=l 

The six patterns are reported at the densities corresponding approximately to the 
following experimental pressures: 0, 3, 5, 7, 10, 15 GPa. In the right-hand panel of 
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Table 1. Parameters used for the form factors [45l l46]. 
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figure [3] we report the calculated neutron diffraction patterns at the same densities; 
these are obtained from the same expressions but with /^(g) replaced with ba, the 
neutron diffraction length [1] . In this figure we also show some experimental data from 
the recent study by Drewitt and co-workers [12] at densities close to the simulation ones. 
These data were collected at the following pressures: 2.2, 4.9 and 8.0 GPa. 

By comparing figure [3] with Fig. 1 in ref. |Tl], it can be readily appreciated that the 
agreement with the experimental X-ray structure factors is excellent at all pressures. 
Our data show that the First Sharp Diffraction Peak position (FSDP) shifts toward 
higher wavenumbers upon pressure increase until it merges with the principal peak, in 
agreement with the data in ref. [11]. From the fact that our simulations are in good 
agreement with the x-ray structure factors from [Tl] as well as the the more recent 
neutron data [12] and the bond-lengths from the most recent EXAFS [8] we conclude 
that the different sets of experimental data themselves are mutually consistent (see also 
discussion below). The difference in the interpretation of the X-ray data by Guthrie et 
al. [H] (i.e. sharp vs. smooth tetrahedral to octahedral structure) is therefore likely to 
be due to the difficulty of extracting good coordination numbers from noisy data. 

The comparison with the neutron data from Guthrie et al. [11] is inhibited by the 
noisiness of the experimental data but the data of Drewitt et al [12] is of higher quality. 
For densities above p ~ 5.0 g/cm^, corresponding to an experimental pressure of 8-9 
GPa, there are, to the best of the authors' knowledge, no available neutron data. Our 
simulations are consistent with the main trends observed in the neutron work, namely 
a diminution in the height and a shift to higher q of the FSDP with increasing pressure. 
However, the agreement between simulations and experiment is less good than in the 
X-ray case. The simulations appear to overestimate the height of the FSDP, especially 
at low pressure, and also to underestimate the amplitude and sharpness of the principal 
peak at about 2.5 A"^. In the total neutron structure factor of Ge02 there is a near- 
total cancellation between the different partial structure factors in the vicinity of the 
principal peak, so that the "principal" peak in the total S{q) is of much lower amplitude 
than in the individual Ge or O partials. The discrepancy we noted above presumably 
reflects a small error in position or width of the principal peak in these partials. The 
FSDP on the other hand reflects the intermediate range order in the glass and the im- 
plication of a disagreement in this domain is that the quenching method we have used 
to prepare the glass does not allow this structure to develop in the same way as in the 
experimental material. In support of this contention, we note that the S{q)s calculated 
from the cold-compressed and densified samples differed most in this region of g-space. 
The position of the FSDP obtained from the neutron structure factors is shown in the 
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Figure 3. Top: simulated X-ray and neutron diffraction patterns as a function of 
increasing density. The experimental data in the right-hand panel (line with error bar) 
are from [12] . Bottom: pressure dependence of the FSDP (extracted from the neutron 
patterns) position obtained from the X-ray diffaction pattern. The experimental data 
(filled green squares) are from |12) . 
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bottom part of figure [3] and increases linearly with density. The corresponding values 
from the Drewitt et al work are also shown. 

Overall then, our simulations agree well with the available pressure- dependent 
diffraction data when we compare directly with the experimentally measured quantity, 
S{q). There is also good agreement with the density dependence of the nearest- 
neighbour separation obtained from analysis of the EXAFS data. The examination of 
the local structure in the simulations supports the conclusion that the transition from 
a tetrahedral to an octahedral glass is smooth and gradual and probably not completed 
even at pressures as high as 15 GPa. 

5. Liquid germania 

The high-pressure behavior of liquid germania presents some differences compared with 
the glassy case. Although a similar increase in mean coordination number is observed, 
both the Ge-0 coordination number and Ge-0 distances (not shown) increase more 
smoothly as pressure is raised. The increase in coordination number begins as soon as 
the pressure in increased. In figure H] we show the oxygen and germanium self diffusion 
coefficients as a function of increasing pressure for liquid germania at 4000 K. It can 
be seen that the diffusion coefficient exhibits the same anomalous maximum in the P 
= 15-25 GPa range as observed in other tetrahedral-network materials (like silica and 
water |5]). In the case of liquid silicates, it was unambiguously shown by Angell et 
al. that such a pressure enhancement of ion mobilities was linked to the formation of 
fivefold coordinated silicon ions, which act as defects in the original tetrahedral network 
structure. [34]. The proportion of GeOs polyhedra with each coordination number is 
shown on figure HI and we can see that its variation with pressure matches exactly the 
variation of the diffusion coefficients, in agreement with the silica situation. 

6. Conclusions 

In this work we used a new polarizable interaction potential to study the high pressure 
behaviour of glassy and liquid germania. For glassy germania, we were able to reproduce 
all the experimental structural information by making comparisons at the same density. 
The density-evolution of the Ge-0 coordination number, bond distance. X-ray and neu- 
tron structure factors was in good agreement with the data from refs. [H [TTl [12]. The 
only observed shortcoming was the inability of our simulations to properly reproduce 
the equation of state for glassy germania. This is most likely due to the relatively short 
time-scale available to computer simulations. From our data and from a comparison 
with the available experimental evidence, it can be concluded that the transition from a 
tetrahedral to an octahedral glass is smooth and gradual and probably not completed at 
pressures as high as 15 GPa. In view of our results, it seems likely that the differences 
between the mechanisms which have been proposed to interpret this data (such as sharp 
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Figure 4. Simulated diffusion coefficients on liquid (T = 4000 K) germania and 
percentage of GeO^ units as a function of increasing pressure. 



vs. smooth tetrahedral to octahedral transitions) is due to the difficulty of extracting 
reliable coordination number from noisy data. Finally, a study of the percentage con- 
centration of Ge04, GeOs, GeOg units shows that as density is increased, these are 
present in different proportions but a state with GeOs units only was not observed. 
As for liquid germania, our results are mainly predictive due to the lack of experimental 
data, and of interest for the comparison with the behaviour of other systems which are 
tetrahedrally coordinated at ambient pressure. We find that liquid germania undergoes 
a transition from a fourfold- to a sixfold-coordinated phase and that this transition is 
smoother than the one observed in the glassy phase. Most importantly, we find that 
the diffusivity in germania behaves anomalously as a function of pressure, showing a 
maximum in the P = 15-25 GPa range, as found in other tetrahcdral-network liquid 
such as silica and water. The cause of this behaviour is traced back to the presence of 
pentahedrally-coordinated Ge ions. 
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